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A computational analysis of a Gulfstream isentropic external compression supersonic inlet coupled to a 
Rolls-Royce fan was completed. The inlet was designed for a small, low sonic boom supersonic vehicle 
with a design cruise condition of M = 1 .6 at 45,000 feet. The inlet design included an annular bypass duct 
that routed flow subsonically around an engine-mounted gearbox and diverted flow with high shock 
losses away from the fan tip. Two Reynolds-averaged Navier-Stokes codes were used for the analysis: an 
axisymmetric code called AVCS for the inlet and a 3-D code called SWIFT for the fan. The codes were 
coupled at a mixing plane boundary using a separate code for data exchange. The codes were used to 
determine the performance of the inlet / fan system at the design point and to predict the performance and 
operability of the system over the flight profile. At the design point the core inlet had a recovery of 96 
percent, and the fan operated near its peak efficiency and pressure ratio. A large hub radial distortion 
generated in the inlet was not eliminated by the fan and could pose a challenge for subsequent booster 
stages. The system operated stably at all points along the flight profile. Reduced stall margin was seen at 
low altitude and Mach number where flow separated on the interior lips of the cowl and bypass ducts. 
The coupled analysis gave consistent solutions at all points on the flight profile that would be difficult or 
impossible to predict by analysis of isolated components. 
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Objectives 

1 . Develop CFD tools for coupled inlet / fan interaction. 

2. Apply the tools to a Gulfstream supersonic inlet coupled to the fan 
for a Rolls-Royce engine. 

3. Determine the performance of the system at the design point. 

4. Assess the operability of the system over the flight profile. 
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AVCS code 


Axisymmetric Navier-Stokes duct analysis code by D. Tweedt and R. Chima 

• Node centered finite-difference formulation 
•AUSM + upwind differencing 

• Axisymmetric blockage term used to model gearbox blockage 

Explicit Runge-Kutta solver 

• Variable At 

• Implicit residual smoothing 
Turbulence model 

•Wilcox 2006 k-co turbulence model and stress limiter 
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SWIFT code 


3-D multiblock Navier-Stokes turbomachinery analysis code by R. Chima 

• Node centered finite-difference formulation 
•AUSM + upwind differencing 

Explicit Runge-Kutta solver 

• Variable At 

• Implicit residual smoothing 
Turbulence model 

•Wilcox 2006 k-co turbulence model with stress limiter 

Coupled to AVCS with a steady mixing plane model 

•SWIFT blade-to-blade solution perturbed about AVCS solution using Giles’ 
characteristic BC 
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Code Cou pling - SYNCEX 




SYNCEX - C language code written by Dr. Dan Tweedt of AP Solutions, Inc. 

• Runs in the background and handles data communication, storage, and 
synchronization between CFD codes 

• User routines read and write boundary condition data to SYNCEX 

• General interpolation routines provided 
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Applicat ion - Gulfstream Quiet 


Supersonic 



•Conceptual design by Gulfstream Aerospace Corporation (GAC) 

• Mach 1 .6 cruise at 45,000 ft 

• Extensive shaping to minimize sonic boom for flight over land 
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Gulfstream Inlet / Rolls-Royce fan 



Gearbox Fairing 


rnal Compression lrtie| 


Flow-Through Struts 


Nozzle Shroud 


Bypass Exhaust 
► Core Exhaust 


Rolls-Royce Engine 


Nozzle 


Bypass Flow 

Core Flow 


Mach 1 


o.o 

Mach No. 


• Isentropic external compression inlet 

• Rolls-Royce engine 

• Plug nozzle 

• Novel bypass duct routes flow internally around large gearbox 
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Computat ional grids 
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SIGG grid code by D. Tweedt 

•Algebraic grid with elliptic smoother used for inlet, bypass, and external flow 
• y+ = 1 _ 2.5 


Region 

Type 

Size (x, r) 

Total points 

inlet 

H 

419 x 95 

39,805 

bypass 

H 

439 x 55 

24,145 

external flow 

H 

360x128 

46,080 

Total 



110,030 
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Gearbox f airing blockage 


™ o.o 

% 0 blockage 




Gearbox fairing modeled using specified blockage in bypass duct grid 
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3-D fan g rid 

TCGRID grid code by R. Chima 
•Algebraic H-grid upstream 

• Elliptic C-grids around blades 
•Algebraic O-grids in tip clearance region 

• y+ = 1 _ 2.5 


Region 

Type 

Size (x, q, r) 

Total points 

upstream 

H 

45 x 30 x 95 

128,250 

rotor 

C 

257 x 46 x 95 

1,123,090 

rotor tip 

0 

201 x 15 x 15 

45,225 

Total 



1,296,565 
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Grid details 



1 . Inlet centerbody radius matched to rotor hub. 

2. Tip clearance is unknown. Used 0.025 inches. 
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Grid details 



Blade to blade grid at tip 
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Inlet Corrected Mass Flow [Ib/sec] 


Solution details 



•Solutions run 10,000 - 20,000 iterations at CFL = 2.5 

• Fan pressure ratio and mass flow monitored for convergence 

• 7 hr. / 10,000 iterations on dual core Xeon CPUs at 3.8 GHz 
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Static pres sure contours 



•Shocks from fan run upstream into diffuser 
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Mach contours 
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Core Recovery Bypass Recovery 


Bypass C apture Ratio and Core Recovery 
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Inlet Capture Ratio 

• Design point: M = 1 .7 (over wing) at 45,000 ft. 

•Spillage = 1.6 % 

• Bypass recovery = 0.885 

• Core recovery = 0.96 
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Fan Pressure Ratio Adiabatic Efficiency 


Fan Effic iency and Pressure Ratio 



390 400 41 0 420 430 440 



• Hub distortion reduces maximum mass flow by 2.4 % 
and peak efficiency by 0.8 points 
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Spanwise profiles at AIP 



•Minimal tip distortion due to bypass design 
•Large hub distortion - uncommon for fans 
• Ps « constant, good BC for isolated inlet calculations 
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Spanwise profiles at fan exit 
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Total Pressure Ratio Adiabatic Efficiency 


• Hub distortion persists through the fan. 
Effects on core compressor are unknown. 
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Nominal Flight Profile 



Mach 

Coupled inlet / fan calculations made at each point shown on flight profile 
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Fan Corrected Conditions 



Mach 
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Boundary Conditions 

Upstream 

•Subsonic: P 0 , T 0 specified, /T extrapolated from interior 
•Supersonic: all quantities specified 

Downstream and bypass exit 

• Bypass closed for M < 0.3 

• Subsonic: Ps specified 
•Supersonic: all quantities extrapolated 

Freestream 

• Characteristic perturbation about freestream 

• Reynolds number varies with altitude: Re = /(P 0oo ,r 0 J 

Fan 

• Corrected speed varies with altitude: Nc = O.Jq 

• Q = constant over most of flight profile 

• Q reduced for M > 1 .4 due to T3 limit based on cycle deck data from GAC 

• Ps = linear with Nc, matched to desired mass flow at 2 operating points 
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Mach Contours 


The following slides show Mach contours at most points on the flight profile. 
• Note boundary layers at fan face, bypass choking, and shock structure. 
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Mach Contours 


M = 0.01, A = 0 FT 



Ground idle 

• Bypass closed 

• Moderate tip distortion 
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Mach Contours 



Takeoff 

• Bypass closed 

• Severe tip distortion 
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Mach Contours 



M = 0.3, 5,000 ft. 

• Bypass open 

• Severe tip distortion 
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Mach Contours 



M = 0.45, 10,000 ft. 

• Diffuser choked 

• Moderate tip distortion 
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Mach Contours 


M = 0.72, A = 18,000 FT 



M = 0.72, 18,000 ft. 

• Bypass nearly choked 

• Low tip distortion 
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Mach Contours 



M = 0.90, 30,000 ft. 

• Bypass choked 

• Shock on shoulder 

• Low tip distortion 
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Mach Contours 


00, A = 30,000 FT 




M = 1 .00, 30,000 ft. 

• Possible shock on nacelle 

• Shock on shoulder 
Low tip distortion 
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Mach Contours 



M = 1 .10, 30,000 ft. 

• Shocks on spike, cowl, shoulder 

• Low tip distortion 
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Mach Contours 
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Mach Contours 



Cruise, M = 1 .60, 45,000 ft. 

• Shocks on spike, cowl, shoulder 

• Hub distortion increasing 
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Mach Contours 
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Capture S tream I i nes 

The following slides show capture streamlines at some points on the flight profile. 

• Note the streamline shape around the bypass splitter lip and the effect on the 
core tip boundary layer. 
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Capture S tream I i nes 
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Capture S tream I i nes 
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Capture S tream I i nes 


M = 0.30, A = 5,000 FT 




M = 0.3, 5,000 ft. 

• Bypass open 

• High negative incidence on bypass lip 

• Streamlines separate at bypass lip 

• Severe tip distortion 
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Capture S tream I i nes 
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Capture S treamlines 



Over wing, M = 1 .70, 45,000 ft. 
• High hub distortion 
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Inlet Core Recovery 


Inlet Core Recovery 

1.00 

0.99 


0.98 
0.97 
0.96 

0.95 
0.94 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 

Mach 

• Inlet recovery is low at M < 0.7 due to bypass separation 

• Recovery is excellent for 0.7 < M < 1 .4 

• Recovery decreases with M > 1 .4 due to shock losses 
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Fan Pressure Ratio 


Fan Pressure Ratio 


2.001 1 1 1 1 1 1 1 r 
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Mach 

• Fan pressure ratio varies with corrected speed 

• Decreases for M > 1 .4 due to reduced Q 
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Fan Pressure Ratio Adiabatic Efficiency 


Fan Pressure Ratio and Efficiency Maps 



380 390 400 41 0 420 430 440 450 



•3 speed lines shown for fan with cruise 
distortion profile. 

• Fan operating line is nearly linear except for 
low Mach numbers. 

•Nc varies between 92 and 107 % of design. 

•M = 0.15 - 0.3 is close to the stall line due to 
high tip radial distortion. 

• Fan PR drops with Nc for M > 1 .4. 

• Fan operates near peak r| at cruise. 
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Conclusions 


Two coupled CFD codes were used to model the interaction between a 
supersonic Gulfstream inlet and a Rolls-Royce fan. 

Advantages of coupled analysis 

• Eliminates modeling at interface 

• Can track inlet distortion effects all the way through the fan 

• Hard to model isolated components for some cases, e.g. low Mach nos. 

Disadvantages 

• Increased grid / solver complexity 

• Longer solutions times 

Gulfstream inlet / Rolls-Royce fan 

• Design point: good inlet recovery, fan operates near peak efficiency. 

• Hub radial distortion passes through fan, may cause problems for the core. 

• Flight profile: stable operation at all operating points. 

• Identified reduced stall margin at low Mach nos. due to bypass lip separation. 
Lip redesign is being studied. 

•2 papers at AIAA 48 th Aerospace Sciences Meeting, Jan. 4-7 2010, Orlando, FL 
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